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Abstract 

We calculate purely weak virtual one-loop corrections to the production cross section of top- 
antitop pairs at the Large Hadron Collider via the gluon-gluon fusion subprocess. We find 
^ ■ very small negative corrections to the total cross section, of order —0.6%, but significantly 

larger effects to the differential one, particularly in the transverse momentum distribution, of 
^ ■ order —5% to —10% (in observable regions). In case of parity-conserving spin-asymmetries 

of the final state, c^aw corrections are typically of a few negative percent, with the exception 
of positive and negative peaks at +12% and —5%, respectively (near where the tree-level 
predictions change sign), while those arising in parity- violating asymmetries (which are 
Ph| identically zero in QCD) are typically at a level of a few permille. 

Oh! 
CD . 

: 1 Introduction 

> ■ 

^ ■ Top quark physics may well be the only context where both accurate Standard Model (SM) 

c3 ' tests and searches for new physics Beyond the Standard Model (BSM) will be carried out 

at the Large Hadron Collider (LHC). If no BSM physics exists at the TeV scale or the 
typical mass scale of new particles is (just) above the energy reach of the machine, one 
may well conceive that most of the experimental and theoretical efforts will concentrate 
in establishing the true nature of the top quark, which in turn will also enable one to 
constrain possible manifestations of new physics. While top quarks have been discovered 
and studied at the Tevatron, the reduced number of events available there will only allow 
one for a percent level determination of the top mass (currently, m t = 172.7 ± 2.9 GeV). 
This precision will be improved by over a factor of two at the CERN machine. Here, one 
will also be able to measure the top-quark width and quantum numbers (i.e., the electric 
charge and isospin, accessible through its Electro- Weak (EW) couplings). While there is 
certainly scope to investigate the EW couplings of top-quarks by resorting to events with 
radiated photons and Z bosons [1] , the V — A structure (or otherwise) of the charged decay 
current can already be probed directly in ti events, if one recalls that the top- (anti) quark 
decays into a bottom- (anti) quark and a W boson rather than hadronising. Finally, for the 
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same reason, the top- (anti) quark transmits its spin properties to the decay products rather 
efficiently, so that the latter can be explored in suitable experimental observables [2, 3, 4]. 

Clearly, in order to perform all the relevant measurements in ti events, any source of SM 
corrections should be well under control. While complete one- loop results exist for QCD [5], 
similar weak effects have been unavailable until very recently [6, 7]. These last two papers 
were concerned with purely weak a\a^ effects entering the qq — > ti subprocess only. It is 
the purpose of this paper to complement those studies, by computing the corrections to the 
gg — > ti channel, which is in fact dominant at the LHC (whereas the quark initiated one 
is the leading partonic component at the Tevatron). Early, though incomplete results, for 
«g«w corrections to top-antitop hadro-production can be found in Ref. [8] for both gg and 
qq initiated subprocesses. Concerning the gg — > ti case, unlike Ref. [8], notice that we have 
also included here the one-loop triangle contributions for gg — > Z* — > ti, which are in fact 
non-zero for off-shell Z bosons. (Recall that the Landau- Yang's theorem [9] is only valid 
for on-shell Z bosons and we have explicitly verified this to be the case in our calculation if 
we take the appropriate limit.) We have also updated the analyses of [8] to the most recent 
top and Higgs mass values as well as Parton Distribution Functions (PDFs). 

The paper is organised as follows. The next section illustrates the importance that EW 
effects should have at TeV energy scales. Sections 3 and 4 will be devoted to describe our 
computation and present the numerical results, respectively. The last section contains our 
conclusions. 

2 EW effects at TeV scale energies 

The purely weak (W) component of Next-to-Leading Order (NLO) EW effects produces 
corrections of the type aw log 2 (/i 2 /M^), where aw = «em/ sin 2 6^, with oem the Electro- 
Magnetic (EM) coupling constant and 6w the Weinberg angle. Here, \x represents some typ- 
ical energy scale affecting the top-antitop process in a given observable, e.g., the transverse 
momentum of either the top (anti)quark or the top-antitop invariant mass. For large enough 
fi values, such EW effects may be competitive not only with Next-to-NLO (NNLO) (as 
a w ~ al) but also wi th NLO QCD corrections (e.g., for \i = 0.5 TeV, log 2 (/i 2 /M^) « 10). 

These 'double logs' are of Sudakov origin and are due to a lack of cancellation between 
virtual and real py-emission in higher order contributions. This is in turn a consequence of 
the violation of the Bloch-Nordsieck theorem in non-Abelian theories [11, 12]. The problem 
is in principle present also in QCD. In practice, however, it has no observable consequences, 
because of the final averaging of the colour degrees of freedom of partons, forced by con- 
finement into colourless hadrons. This does not occur in the EW case, where the initial 
state generally has a non-Abelian charge, as in proton-proton scattering. Besides, these 
logarithmic corrections are finite (unlike in QCD), since Myy provides a physical cut-off for 
iy-emission. Hence, for typical experimental resolutions, softly and collinearly emitted weak 
bosons need not be included in the production cross section and one can restrict oneself to 
the calculation of weak effects originating from virtual corrections only. By doing so, similar 
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logarithmic effects, ~ aw log 2 (/i 2 /M|), are generated also by Z-boson corrections. Finally, 
in some instances all these purely weak contributions can be isolated in a gauge-invariant 
manner from EM effects which therefore may not be included in the calculation (as it is the 
case here). (Besides, EM corrections are not subject to Sudakov enhancement.) In view 
of all this, it becomes of crucial importance to assess quantitatively such weak corrections 
affecting, in particular, a key process (for both present and future hadron colliders) such as 
top-antitop hadro-production. 



3 Calculation 



It is the aim of our paper to report on the computation of the full one-loop weak ef- 
fects entering the subprocess gg — > ti, through the perturbative order agOw- We will 
instead ignore altogether the contributions of tree-level a\a^ terms involving the radiation 
(bremsstrahlung) of real Z bosons. In the Feynman-'t Hooft gauge 1 , the one used for this 
calculation, neglecting the 6-mass, one has to calculate the following one-loop prototype 
diagrams for gg — * ti, ignoring permutations of external gluons: 



' '-^-i-rr&Wi— 






Here, we mean that the numerator of the massive gauge boson propagators is taken to be —ig^u and 
Goldstone bosons, with masses equal to their gauge boson counterparts, are included where appropriate. 
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(In the calculation of all such graphs, we have retained the full spin dependence of the 
final state particles.) 

Given the large number of diagrams involved in the computation, it is of paramount 
importance to perform careful checks. In this respect, we should mention that our expres- 
sions have been calculated independently by at least two of us using FORM [13] and that 
some results have also been reproduced by another program based on FeynCalc [14]. Upon 
removing the one-loop triangle contributions for gg — > Z* — > ti and for identical choices of 
Higgs mass, we also reproduce well the results of the first paper in Ref. [8]. (In fact, this 
triangle contribution is always only marginally relevant.) Finally, we find reasonable agree- 
ment with Ref. [4] (see also [15, 16]) in the Sudakov limit, i.e., for large invariant masses 
and transverse momenta of the final state, provided that the final state particles are fairly 
central (see later on). 

Some of the diagrams contain ultraviolet divergences. In the case of self-energies for 
massive particles the mass subtraction has been effected on mass-shell so that the masses 
refer to the physical (pole) masses. The remaining divergences have been subtracted using 
the 'modified' Dimensional Reduction (DR) scheme at the scale fx = Mz- The use of DR, as 
opposed to the more usual 'modified' Minimal Subtraction (MS) scheme, is forced upon us 
by the fact that the W and Z bosons contain axial couplings which cannot be consistently 
treated in ordinary dimensional regularisation. The strong coupling is not renormalised 
by the weak interactions, which means that there are Ward identities which cancel the 
divergent corrections to the strong coupling. Thus the choice of subtraction scheme has no 
effect on our final results since the scheme dependence cancels when all graphs are summed 
over. On the other hand the the EM coupling, oem, has been taken to be 1/128 in order to 
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correctly account for the SM running of the electroweak coupling up to the threshold for ti 
production. 

For the top mass and width, the latter entering some of the loop diagrams, we have 
taken m t = 175 GeV and T t = 1.55 GeV, respectively. (As already intimated, the 6-quark 
was considered massless.) The Z mass used was Mz = 91.19 GeV and was related to the W 
mass, Mw, via the SM formula Mw = Mz cos 9w, where sin 2 9w = 0.232. (Corresponding 
widths were Y z = 2.5 GeV and IV = 2.08 GeV.) The Higgs boson mass and width were set 
to 150 GeV and 16 MeV by default, respectively. However, other mass (and consequently 
width) choices (above the LEP limit of Mh -115 GeV) have been investigated (see later 
on). The PDFs we have used are CTEQ6L1 [17] taken at the factorisation scale Q = 2mt. 
We have also checked other sets, but found no significant difference in the relative size of 
our corrections. 



4 Numerical results 

Our initial findings are presented in Figs. 1-2. Here, we consider both differential spectra of 
some kinematic observables as well as global asymmetries, the latter plotted, e.g., against 
the invariant mass of the ti pair. Notice that this last quantity can only be defined when 
both the top and anti-top four-momenta are reconstructed, which happens in the case of 
fully hadronic and semi-leptonic/hadronic decays 2 , but not for fully leptonic ones (where 
two neutrinos escape detection). 

The definitions of the asymmetries are as follows: 

AiLdo = da ++ — da + _ + da — da |_, 

A L da = da^ — da + , 
Apv da = da — da ++ . (1) 

For Al only the polarisation of either the t-quark or t-antiquark is assumed to be measured, 
whereas the other two asymmetries require the determination of the polarisations of both 
the outgoing particles. All is parity-conserving while the other two are parity- violating 3 . 
Here, the indices + and — refer to the helicities of right (R) and left (L) handed (anti)top 
quark, respectively. (Other basis choices are also possible, see Ref. [2].) 

We find that the overall effect of our C(ag«w) corrections is about —0.6% at the inclusive 
level (i.e., to the total cross section, as obtained from the integral of any of the curves in 
Fig. 1). However, for differential cross sections, effects can be of either sign, notably in the 
(anti)top transverse momentum and top-antitop invariant mass. For an LHC integrated 
luminosity of 300 fb" 1 , differential rates of order a few 10~ 5 pb may yield detectable events 
(after accounting for decay fractions, tagging efficiency and reconstruction performance). 

2 In the second case one would reconstruct the longitudinal neutrino momentum by equating the trans- 
verse one to the missing transverse energy and enforcing the W mass reconstruction. 
3 See Ref. [2] for a choice of observables correlated to these asymmetries. 
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In the corresponding observable kinematic range, the maximum correction occurs for the 
transverse momentum spectrum of the (anti)top at around 1 TeV (in the Sudakov limit), 
where it reaches almost the —10% level. In the same kinematic regime, the effects are smaller 
for the invariant mass and pseudorapidity distributions. At small transverse momentum, 
invariant mass as well as in the very forward/backward direction is where the corrections 
are positive, with a maximum of 0(+2%) for the second of these observables. 

In the case of the parity-conserving asymmetry, for which the 0{a\) result is non-zero, 
O(agttw) effects enter significantly (up to the +12 and —5% level or so) only near the point 
where tree-level predictions are zero. Otherwise, they amount to a few negative percent at 
the most. For the parity- violating asymmetries, the relevant quantity is the actual value of 
the (9(agaw) result, as the O(o§) term is identically zero. In both cases, the rates are at 
the permille level (away from the M t i ~ 2mt threshold, where our fixed order results are 
not fully reliable). 

The dependence on the actual value of the Higgs mass is generally negligible at both 
inclusive as well as differential level, with the possible exception of the pseudorapidity dis- 
tribution in the very central region, see Fig. 3 (top-left frame), where the absolute size of 
the C(agaw) corrections is shown for Mh = 150 and 200 GeV. In fact, the inclusive correc- 
tion varies from —2.33 to —2.51 pb, respectively, in comparison to a (Higgs independent) 
tree-level result of 384 pb. For the other differential spectra studied such effects are rather 
uniformly spread across the given kinematic range. As for the asymmetries, here the effect 
of an increased Higgs mass varies significantly with M t t, yielding differences with respect to 
the rates obtained with our default Higgs mass value which are not negligible over most of 
the kinematical intervals considered, see Fig. 3 (top- right and bottom frames). 

Before closing, it is of interest to compare our exact 0(ag«w) results with those of 
Ref . [4] . Notice that the latter only include the case of opposite helicities 4 in the final state 
and are limited to the contribution of non-angular single (~ «w log(s/M^)) and double 
(~ aw log 2 (s/M^/)) logarithms. In performing the comparison between our results and 
those in Ref. [4] we have removed the EM component from the latter. From Fig. 4, it is 
clear that for the logarithmic approximation described be valid all Mandelstam variables 
s, t, u must be very large, condition which is obviously not fulfilled at small/large scattering 
angles. However, it should be appreciated that the ratio between the full 0(a^a^) term and 
the one obtained in NLO Sudakov approximation is a constant to a very good approximation 
already at moderate energies and for most angles. Thus, in principle, a parameterisation of 
this constant as a function of the angle should be possible, in view of high statistics Monte 
Carlo simulations. 



The contribution due to identical final state helicities becomes not negligible near the threshold at 
M a « 2m t . 
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Figure 1: Differential distributions of the subprocess gg — *> ti through the C?(a§) (top frames, 
dotted) and the C(ag«w) (top frames, solid) as well as the percentage of the latter with 
respect to the former (bottom frames, solid) for the (anti)top transverse momentum p?, the 
top-antitop invariant mass M t i and the (anti)top pseudorapidity r\ t . (Lightly/Red coloured 
solid tracts in logarithmic scale are intended to be negative.) 



5 Conclusions and outlook 



For kinematic variables accessible at the LHC, purely weak corrections through one-loop 
level (without Z bremsstrahlung) to the top-antitop cross section via gluon-gluon fusion 
are generally small, although - in order to obtain both the appropriate normalisation and 
shape of the theoretical prediction - they cannot be neglected in the Sudakov regime of 
some observables. In contrast, in line with the results reported in [18]-[21] for the case of 
massless quark pair production, such one-loop weak effects are always crucial in massive 
quark pair production when spin-asymmetries are considered (particularly, parity-violating 
ones). Ultimately, our results will have to be put together with those from Refs. [6, 7] 
(which we are in the process of repeating), in order to study the top-antitop cross section 



7 




2000 3000 



3000 3000 
M„ (GeV) 



3000 3000 
M„ (GeV) 



Figure 2: The differential spin asymmetry All (as defined in the text) of the subprocess 
gg — > tt through the 0(a\) (top frame, dotted) and the C(ag«w) (top frame, solid). (Note 
that the LO QCD contribution changes sign at « 900 GeV and is heavily dependent on M t i 
whereas the 0(a^a^) correction is not.) Just below the top frame we show the percentage 
correction to the (non-zero) LO QCD asymmetry for A L l due to 0{a\a^j) effects. The 
lower two frames display the asymmetries Al and Apy (as defined in the text), which 
vanish exactly in LO QCD, through the same order. The asymmetries are calculated along 
the helicity function of the top-antitop invariant mass M t i. 



at the level of precision required by the LHC experiments. However, particular care should 
eventually be devoted to the treatment of real Z production and decay in the definition of 
the inclusive data sample, as this will determine whether (possibly positive) tree-level Z 
bremsstrahlung effects have to be included in the theoretical predictions through 0{a\a^\ 
which might counterbalance the negative effects due to the one-loop Z exchange estimated 
here. Along the same lines, it should be recalled that NNLO terms ought to be investigated 
too, as it is well known from the Sudakov treatment that they may well be sizable in 
comparison to the NLO ones (see, e.g., Ref. [22]). Finally, we have verified that the effects 
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Figure 3: The absolute size of the 0(a\ass[) corrections to the subprocess gg — > ti for 
the distribution in (anti)top pseudorapidity t] t (top-left frame) and the differential spin 
asymmetries (as defined in the text), for M H = 150 GeV (solid) and M H = 200 GeV 
(dotted). The asymmetries are calculated along the helicity function of the top- 

antitop invariant mass M#. 



studied here for the gg — > ti channel are of no phenomenological relevance at the Tevatron. 
The experimental impact of 0(a^a^) effects will eventually have to be assessed in a proper 
detector simulation, in presence of top-antitop decay, parton shower and hadronisation: 
in fact, as shown in [23], the possibility of extracting such effects is generally limited by 
systematics rather than statistics. 
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Figure 4: A comparison of the exact (9(ag«w) corrections to those obtained from angular 
independent (double and single) logarithms only in Ref. [4], for the two combinations with 
opposite helicities in the final state. The graphs on the left (right) represent the case of 
large (small) angle scattering. Notice that we have subtracted the EM contributions from 
the formulae in Ref. [4]. 
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Abstract 

This is an Erratum to a paper of ours, Phys. Lett. B 639 (2006) 513. After its publication, 
we have discovered a mistake in a numerical program that affects the results presented 
therein. We provide here the corrected version. 

New numerical results 

The numerical program used to produce the results of Ref. [1] was affected by a mistake, 
which has now been corrected. We present in Figs. 1-3 the amended results, in correspon- 
dence to the same figures in the original paper. We also note, as remarked in [2], that the 
parity-violating asymmetry A PV defined in [1] is identically zero through (9(aga w ), as we 
could now verify explicitely. The inclusive results are now as follow: while the (Higgs mass 
independent) tree-level cross section for gg — > ti at the LHC is 384 pb, we now find that 
the weak corrections through 0(a^a^) amount to —9.65 pb for Mh = 150 GeV and —9.39 
pb for Mh = 200 GeV. All these plots and cross sections correspond to the numerical setup 
(masses, couplings, PDFs, scales, etc.) declared in [1]. 

For all such results, as well additional ones presented in Refs. [2, 3] (where independent 
calculations of the same corrections tackled in [1] were performed), we found very good 
agreement between ourselves and Refs. [2, 3], for the same choice of input parameters. 

As for Fig. 4 of [1] , here the differences are not substantial. Besides, the purpose of those 
plots was to illustrate the difference between the full results and the Sudakov approximation, 
which is largely unaffected by the correction we made to our program. Hence, we do not 
reproduce those plots here. 

Finally, we have now also calculated the 0(a\a-w) corrections to the qq — > ti subprocess 
and the results obtained (not shown here) are in agreement with the corresponding ones in 
Refs. [4, 5]. 
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Figure 5: Differential distributions of the subprocess gg — > tt through the 0(a\) (top frames, 
dotted) and the C(agaw) (top frames, solid) as well as the percentage of the latter with 
respect to the former (bottom frames, solid) for the (anti)top transverse momentum px, the 
top-antitop invariant mass M t i and the (anti)top pseudorapidity r\ t . (Lightly/Red coloured 
solid tracts in logarithmic scale are intended to be negative.) 
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Figure 6: The differential spin asymmetry All (as defined in [1]) of the subprocess gg — > tt 
through the 0{a\) (left plot, top frame, dotted) and the C(agQ;w) (left plot, top frame, 
solid). (Note that the LO QCD contribution changes sign at ps 900 GeV and is heavily 
dependent on M t i whereas the 0{a\a^) correction is not.) In the left plot, bottom frame, 
we show the percentage correction to the (non-zero) LO QCD asymmetry for All due to 
0(a\otw) effects. The right plot displays the asymmetry Al (as defined in [1]), which 
vanishes exactly in LO QCD, through the same order. The asymmetries are calculated 
along the helicity axis as a function of the top-antitop invariant mass M t j. 
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Figure 7: The absolute size of the 0(a^a^) corrections to the subprocess gg — > tt for the 
distribution in (anti)top pseudorapidity r\ t (top plot) and the differential spin asymmetries 
(bottom plots, as defined in [1]), for M H = 150 GeV (solid) and M H = 200 GeV (dotted). 
The asymmetries are calculated along the helicity axis as a function of the top-antitop 
invariant mass M t f. 
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